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We study the creation of solitons from particles, using the X(f> 4 model as a prototype. We consider 
the scattering of small, identical, wave pulses, that are equivalent to a sequence of particles, and 
find that kink-antikink pairs are created for a large region in parameter space. We also find that 
scattering at low velocities is favorable for creating solitons that have large energy compared to the 
mass of a particle. 



A wide variety of systems, ranging from polyacetylene 
and Josephson junctions to high energy particle physics 
models, contain non-perturbative, "soliton" or "solitary 
wave" excitations in addition to pertubative "particle" 
excitations [l[ . An important unsolved problem is to find 
ways to transition from the particle sector to the soliton 
sector. At a pragmatic level, we would like to develop 
implcmcntablc schemes that might enable solitons to be 
built out of particles. A transition from two energetic 
particles to solitons, however, is known to be exponen- 
tially suppressed e.g. (for a review see 0) though it 
may occur more readily in certain situations, such as in 
the background of a pre-existing kink,0, Hj]. 

In this paper, we will determine a class of initial condi- 
tions that consist of small amplitude perturbations that 
scatter and successfully lead to the production of a kink- 
antikink ("kk") pair in 1+1 dimensions. A trivial scheme 
to determine such a set of initial conditions is to time 
reverse the annihilation of kk. Then the time reversed 
particles would assemble into an outgoing kk. However, 
in any practical setting, such initial conditions would be 
very hard to arrange since the characteristics of the radi- 
ation from kk annihilation arc highly non-trivial. Instead 
we want to consider "clean" initial conditions in which 
we scatter identical wave pulses, somewhat like 2 parti- 
cle scattering. The simplicity of the initial state comes 
with a price in that the final state will now not only con- 
tain a kk but also some radiation. Our approach thus 
differs from other studies which generally considered ini- 
tial conditions containing a single kink and hence had 
non- vanishing topological charge, e.g. @. 

In order to determine clean initial conditions that give 
kk in the final state, we draw lessons from the sine- 
Gordon model which contains both particle and soliton 
sectors and has been studied extensively, both classically 
and in quantum theory While the complete in- 

tegrability of the sine-Gordon model permits many exact 
solutions, it also leads to a disappointing disconnection 



between the particle and soliton sectors, not present in 
many other models which admit solitons. For example, 
in the sine-Gordon model, it is not possible to start with, 
say, a soliton and an antisoliton and end up with parti- 
cles. If a soliton and an antisoliton are set up to collide 
and possibly annihilate, they simply pass through each 
other. Thus soliton scattering states do not convert to 
particle states (even in quantum theory). 

What is important for us is that the sine-Gordon model 
also contains "breather states". If a breather state has 
large amplitude, it can be interpreted as a bound state 
of a soliton and an antisoliton, in which the two keep os- 
cillating about each other but never annihilate. On the 
other hand, small quantized breathers have been inter- 
preted as fundamental particles in the theory. Then one 
might expect the breather to be a bridge between the 
particle and soliton sectors. In the sine-Gordon model, 
however, the breather is a stable object in itself and fails 
to connect the particle and soliton sectors. 

To connect the particle and soliton sectors it is nec- 
essary to depart from the sine-Gordon model. The 
smaller the departure, the weaker will be the connec- 
tion between the particle and soliton sectors. Then, 
if we depart weakly from the sine-Gordon model, we 
expect long-lived "breather-like" states that can tran- 
sition to both widely separated kink-antikink pair and 
also to particles. Such long-lived states have been 
discovered in various systems and have been termed 
"bions" in certain contexts and "oscillons" in others 
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Motivated by these considerations, we study kk pro- 
duction in the \cj) model 



where we have rescaled fields and coordinates so that 4>o 
is the only parameter in the model. The equation of 
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motion is 
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where overdots denote time derivatives and primes de- 
note spatial derivatives. The mass of a fundamental ex- 
citation can be found by considering small fluctuations 
around one of the vacua (say <f> = +1) and is m = v2- 
The kink profile is 
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The energy of a kink is found from the energy expression 
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E = 4>1 I dx 
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and is 
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Note that the kink energy may be made very large com- 
pared to the particle mass by taking large values of 4>q. 
However, O itself does not enter the classical dynam- 
ics of the scalar field though it does play a role in the 
quantized model. 

We would like to use breather-like solutions in the A</> 4 
model in our initial condition. However, such solutions 
are not known analytically. Hence we simply use the 
breather solutions of the pure sine-Gordon model 
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which are given by 

(j)f(t,x;uj,v) = 
where 
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T = j[t - v(x - x )\ , X = j[x 
{l-v 2 )- 1 ' 2 



x 



vt] 



1 



J] = v 1 — oj 2 /oj . 
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In Eq. ([7|), the tan _1 (-) function is taken to lie in 
the interval (— 7r/2, +7r/2). Apart from the boost 7 and 
shift xq, a breather solution is labelled by the parameter 
oj G (0, 1). The solution for small oj can be viewed as 
a sine-Gordon kink and an antikink that are oscillating 
back and forth, merging and emerging forever. Note that 
the breather is localized around one vacuum (at —1), and 
probes the second vacuum (at +1) for durations that vary 
inversely with oj. For ui « 1, the breather describes oscil- 
lations in the vacuum around <fi = — 1. In the quantum 
theory, these oscillations are quantized and the energy of 
the lowest quantum state is equal to that of a particle, 
leading to the identification of the lowest energy breather 
with the particle excitation in the model. 



Before proceeding consider an initial unboosted sine- 
Gordon breather (Eq. ([7]) with v = 0), in the A0 4 
model with equation of motion given in Eq. ((2|). (That 
is, the initial condition is (f>(0, x) = <$j G (0, x; oj, 0) and 
<j)(0, x) = <j>l G {0, x; oj, 0).) The energy of the solution can 
be obtained by evaluating Eq. (|4]) at t = when <fib = — 1 
for all x. Then the potential and gradient terms do not 
contribute, and the kinetic contribution is easily evalu- 
ated. As in the sine-Gordon model we find 



E h = ^l-oj^ 2 



The ratio of kink to breather energy is 
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where in the last expression we assume oj « 1. The field 
profile itself, (j)(t, x), can be obtained numerically and we 
have checked that it is oscillatory and long-lived. More 
specifically, we have shown that half of the initial energy 
i?b in the simulation box (itself much larger than the 
breather size) is radiated in a time Ti/ 2 — 5 x 10 4 A -19 , 
independently of oj. 

We now turn to the problem at hand, namely the cre- 
ation of kk from particles. Our initial conditions will 
consist of a train of little {i.e. oj « 1) breathers com- 
ing in from the left and another identical train of Nj, 
breathers coming in from the right. We will study the 
collision of these breather trains for a variety of parame- 
ters and look for the formation of kk. Hence our initial 
condition corresponds to an incoming state 
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with 



j[t- v n (x-x 0n )} , X n = ^[(x-x 0n ) -v n t] (12) 



where xo n = a + nd, v n = —v < for n > 0, and 
XQ n = —a + nd, v n = +v > for n < 0. The parameter 
a is half the separation between the trains at t = and d 
is the separation between different breathers in the same 
train. The initial conditions (at t = 0) are 



(j>{0,x)= f{0,x) , ^(0,x) = f(0,x) 



(13) 



To further motivate our choice of initial conditions, 
let us consider what might be required to form a kk. 
Initially, the field is oscillating about the 4> = — 1 vac- 
uum. To form kk, we need the oscillations to extend 
into the <f> = +1 vacuum. So we need to build up the 
field oscillations. As the leading breathers in the trains 
collide, the field at x = starts oscillating. Subsequent 
breathers provide additional kicks to the oscillations at 
x = 0. Provided the subsequent collisions arc in phase, 
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the amplitude of oscillations at x — will grow in res- 
onance. The growth must compete with the dissipation 
due to the emission of particle radiation. If the growth 
wins, the oscillations will eventually extend up to the 
<f> = +1 vacuum and then it will be energetically favor- 
able for <fi(t, 0) to stay there. Then it becomes likely that 
kk will be created. 

The same heuristic argument may be applied to the 
pure sine-Gordon model and serves to show its lim- 
itations. We know that kk are not created in the 
sine-Gordon model, but it is not because dissipation is 
stronger than resonant growth. Instead the integrability 
ensures that the breather trains pass unscathed through 
each other. So the heuristic argument should be taken 
as motivation but cannot be taken too literally; instead 
we must solve the equations of motion and check for kk 
production. However, what seems clear is that kk pro- 
duction may proceed via a resonance and, just as a child 
can swing higher and higher by timing her movement to 
within a factor of 2 per kick, this level of tuning may be 
all that is needed to produce kk. 

The equation of motion for the scalar field, Eq. @, is 
solved numerically with the initial conditions in Eq. 
using the iterated Crank-Nicholson method with two it- 
erations , and absorbing boundary conditions [ifl at 
the ends of the lattice. The fields are evolved for one 
light crossing time. As an additional check, we have also 
evolved the intial conditions using Mathematica, though 
with fixed boundary conditions (see 2J] for the note- 
book). The Mathematica results are generally consis- 
tent with the Crank-Nicholson method but there are a 
few discrepancies. These may be due to the different in- 
tegration routines which also have different accuracies. 
The Crank-Nicholson implementation is more transpar- 
ent and we find it more reliable, while the Mathematica 
implementation is more convenient to use. 

The problem contains many parameters, all related to 
the choice of initial condition: oj, v, a, d, and Nf,. For 
a given value of ui, Eq. dTUJ) shows that, just on ener- 
getic grounds, we need iVj > 0.6/Vl — uj 2 . We have 
taken Nb = int[2/vl — uj' 2 ] + 1 where int[x] denotes the 
largest integer less than or equal to x. Somewhat arbi- 
trarily, we take the initial half-separation of the trains 
to be a — lO/rju = 10/vl — uj 2 , corresponding to 10 (at 
rest) breather widths. The separation of the breathers 
in a train is taken to be d = 2/ yl — uj 2 . The only pa- 
rameters left to specify are uj and v. For a breather to 
have energy comparable to the particle mass, and the 
kink to have energy much larger than the particle, we 
require uj to be very close to 1. With u = 0.99, the kink 
energy is about 4 times that of a breather. We shall take 
uj G (0.90,0.99). We then do runs for different values of 
v and look for kk formation. 

An example of kk production is shown in Fig. [1] where 
we give two snapshots of the evolution. Animations of 
the evolution may be found in Ref. 24|. Generally, by 
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1: Two snapshots of the collision of breather trains for 
w = 0.99, v — 0.43 and other parameters as described in 
the text. T denotes a light crossing time. The initial state 
contains the train of breathers. Subsequently, kinks appear 
and move apart. 



looking at the field profile, it is quite clear when a kk has 
been created. However, there are some instances in which 
the outcome is not so clear-cut. This includes the case 
when the field profile shows kk that are not separated 
by a large distance or are almost at rest with respect 
to each other. Then there is the possibility that the kk 
will annihilate. In such cases, we have chosen to call it 
a kk creation event if the kinks survived for at least the 
duration of the simulation. Another novel outcome we 
have seen is that for some parameters two or more pairs 
of kk are produced. 

In Fig. [2] we plot the region on the (uj,v) plane for our 
choice of parameters that lead to the formation of a pair 
(or more) of kk. Note the trend - higher uj i. e. smaller 
breather energy, requires lower incoming velocity. This 
indicates that it is preferable to scatter many particles at 
low energy to create solitons. If the incoming velocity is 
too high, the breather trains simply pass through, as in 
the sine-Gordon model. Also, note the occasional holes 
in the plot (e.g. oj = 0.91, v ~ 0.82) where we did not 
observe kk formation. This substructure in the plot is 
reminiscent of the bands observed in kk scattering [22| 
and suggests that kk formation may be due to resonance. 

The region leading to kk formation is reasonably large 
but does not extend to arbitrarily high uj. For example, 
we have not found initial conditions leading to kk for- 
mation for uj > 0.99. The expanse of the "successful" 
region does not concern us at the moment because our 
main objective was to find a set of clean initial conditions 
that led to the formation of kk. We would be surprised 
if future investigations do not find a larger set of suc- 
cessful clean initial conditions, even for very high values 
of uj. Whether these initial conditions are achievable in 
a practical setting is a separate matter, and depends on 
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FIG. 2: Results from our Crank-Nicholson code (square sym- 
bols) mapping out kk formation in the (u, v) plane for the 
choice of other parameters as described in the text. Note the 
occasional gaps where kk are not formed, and the downward 
trend with larger uj (weaker incoming pulses). 



the details of the experiment. 

There are several directions in which it would be useful 
to extend our results. The first is to scan the space of 
initial conditions more carefully, to gain further under- 
standing of what conditions enable kk formation. Our 
space of initial conditions could also be enlarged, if nec- 
essary. For example, different breathers in a train could 
come in with different velocities. We could also envi- 
sion "building up" by starting with very large ui (small 
energy) breathers, and building states corresponding to 
smaller lu (larger energy) , which can then collide to form 
kk. Another direction is to include quantum effects in 
the scattering. This would require more precise under- 
standing of the breather and kink states in terms of par- 
ticles. In the quantum sine-Gordon model, soliton op- 
erators have been written down in terms of an infinite 
number of particle operators (26|. We expect that the 
soliton operator in the \<p 4 model should be expressible 
in terms of a finite number of particle operators oth- 
erwise it would seem impossible to build a kk starting 
with particles. Yet another direction to proceed would 
be to consider solitons in higher dimensions. Then we 
can study the creation of vortex-antivortex or monopolc- 
antimonopole pairs in suitable systems. We would clearly 
need higher dimensional analogs of breathers and we ex- 
pect that oscillon states can play this role. Finally, it 
would be useful to generalize our initial state to real sys- 
tems. After all, polyacetylene is described by the A</> 4 
model and we may expect to be able to create kk there. 
(Similar problems also arise in polymer physics in the 



context of polymers that pass through a membrane [271].) 
Our results do not directly apply to polyacetylene be- 
cause the dynamics there is non-relativistic. However, 
with suitable generalization, it may become possible to 
test some of these ideas experimentally. 
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